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ABSTRACT 



X 



We investigate the effect of pulsar polar cap (PC) heating produced by positrons 



o 
o 

CN ' returning from the upper pair formation front. Our calculations are based on a 

^ • self-consistent treatment of the pair dynamics and the effect of electric field screening 

. by the returning positrons. We calculate the resultant X-ray luminosities, and 

^ . discuss the dependence of the PC heating efficiencies on pulsar parameters, such as 

characteristic spin-down age, spin period, and surface magnetic field strength. In this 

^ , study we concentrate on the regime where the pairs are produced in a magnetic field by 

' curvature photons emitted by accelerating electrons. Our theoretical results are not in 

■ conflict with the available observational X-ray data and suggest that the effect of PC 

Tjj- ■ heating should significantly contribute to the thermal X-ray fluxes from middle-aged 

. and old pulsars. The implications for current and future X-ray observations of pulsars 

Q . are briefly outlined. 

Q ■ Subject headings: pulsars: general — radiation mechanisms: nonthermal — relativity 

> ■ — stars: neutron — X-rays: stars 



1. INTRODUCTION 

X-ray emission has been detected from several dozen pulsars in observations by ROSAT 
(Becker &: Trumper 1997), ASCA (Saito 1998) and Chandra (Zavlin et al. 2000). In cases where 
good spectral measurements are available, the emission in the X-ray band seems to have both 
thermal and non-thermal components. While the non-thermal components are most likely of 
magnetospheric origin, especially when they connect smoothly to higher energy spectra, the origin 
of the thermal components is less clear because there are several mechanisms for production of 
thermal emission in the soft X-ray band. Pulsars younger than about 10^ yr are expected to 
have significant neutron star (NS) cooling components and observed X-rays from these pulsars 
may be consistent with standard cooling scenarios (Page 1998). However, many older pulsars, 
including millisecond (ms) pulsars, beyond the age where cooling emission drops sharply, have 
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strong thermal X-ray emission. In this case, some heating mechanism(s) must be operating. Two 

which have been proposed are heating by particles, accelerated in the magnetosphere, flowing back 
to the NS PCs (PC heating: Ruderman & Sutherland 1975, hereafter RS75; Arons 1981, hereafter 
A81) and internal heating by frictional forces in the crust (Shibazaki & Lamb 1989). 

Because PC heating emission is intimately tied to the pulsar acceleration mechanism. X-ray 
emission of older pulsars provides very strong constraints on pulsar models. Models which predict 
X-ray luminosities from PC heating which are in excess of observed luminosities are not viable. 
The earliest PC models, based on vacuum gaps at the poles (RS75), predicted PC heating 
luminosities that were way above the first Einstein detections and upper limits on pulsars (Helfand 
et al. 1980). However, predictions of the space-charge limited flow (SCLF) model of A81 were 
within all observed values. Similarly, the PC heating predictions of the first outer-gap models 
(Cheng, Ho &: Ruderman 1986), where particles are accelerated in the outer magnetosphere, were 
well above the observed X-ray luminosities. Later models allowed for a reduction of this huge 
PC X-ray emission by scattering in a "pair blanket" which forms near the surface (Wang et al. 
1998). It is clear that models with PC vacuum gaps and outer gaps, where half of all the particle 
acceleration energy returns to the surface to heat the PCs, will have severe problems preventing or 
hiding the excess X-ray emission. In this respect, SCLF models have a distinct advantage in that 
the number of particles that are accelerated downward to heat the PC are only a small fraction 
of the number of primary particles that are accelerated upward. This is because the accelerating 
electric field (parallel to the magnetic field) arises not from a pure corotation charge, but from a 
slight imbalance between the corotation charge and the actual charge. The field grows with height 
above the surface until electron-positron pairs from the accelerated, upward moving electrons 
screen the field above a pair formation front (PFF). A81, in his study of SCLF acceleration, found 
that the screening is primarily due to positrons which turn around and accelerate backward to the 
stellar surface. The number of returning ("trapped") positrons required to screen the is very 
small compared to the number of primary particles (electrons). In his studies, he assumed that 
pairs were produced only by curvature radiation (CR) photons. 

We have begun an updated study of PC SCLF acceleration (Harding & Muslimov 1998; 
hereafter HM98) based on the work of Muslimov & Tsygan (1992; hereafter MT92), who found 
that including effects of general relativity produced a substantial difference in the parallel electric 
field induced above the surface of a NS. In particular, the effect of inertial frame dragging induces 
a field that is typically a factor of 50 times higher than that derived by A81, which is that induced 
by flaring of polar field lines alone. This relativistic effect of inertial frame dragging and the 
consideration of radiation from inverse Compton scattering (ICS) of thermal X-rays from the NS 
surface by the primary electrons introduces profound changes in SCLF models. We (HM98) found 
that pairs are produced by ICS photons from accelerating electrons at lower altitudes than are 
pairs from CR photons. Pairs from ICS radiation could therefore screen the E'y before any CR 
photons are produced. We further found that positrons accelerating downward from ICS PFFs 
could potentially cause unstable acceleration if ICS pairs from positrons screened the E\\ above 
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the surface. In our previous study, we computed self-consistently the height at which screening 
begins by assuming that the CR and ICS PFFs are located at the point where the first pair is 
produced. Wc also have given rough estimates for the returning positron flux and polar cap 
heating luminosity expected from CR PFFs (Zhang & Harding 2000). However, we did not model 
the structure of the PFF and thus could not determine the flux of returning positrons as a fraction 
of the flux of primaries (the returning positron fraction) or even whether full screening occurs. 
Many unanswered questions therefore remain, such as whether ICS photons can always screen the 
£^11 and whether there are enough returning positrons to produce pairs which screen E'n near the 
surface. 

To address these questions, as well as provide PC heating rates for the new SCLF models, we 
have made a detailed and self-consistent study of the screening at both CR and ICS PFFs. We 
have found that while ICS photons arc able to produce PFFs in nearly all pulsars, ICS radiation is 
able to completely screen the i^y only in pulsars with higher surface magnetic fields. The magnetic 
field value required for ICS screening is strongly dependent on the surface temperature of the NS. 
For PC temperatures T < 10^ K, pulsars with surface Bq ~ 0.1 Bcr, where B^r = 4.4 x 10^^ G 
is the critical field strength, are capable of ICS screening, but for T ~ 10^ K, the required field 
decreases with temperature. If the ICS photons are unable to completely screen the E'n , then the 
primary electrons will keep accelerating until they radiate CR photons which can produce a PFF. 
In some cases, such as ms pulsars, incomplete ICS screening via non-resonant scattering can make 
a significant contribution to PC heating, especially for those pulsars where CR photons cannot 
produce any pairs. The results of these studies will be presented in two papers. This first paper 
will discuss CR screening and PC heating in lower field pulsars, i.e. those below the threshold 
for ICS screening. In a second paper, we will discuss ICS screening in higher field pulsars, as well 
as partial ICS screening in older and ms pulsars, and will address the issues of lower PFFs and 
stability of ICS PFFs. In the present paper, we first discuss the SCLF electric field solution used 
and the conditions for the screening of this field. We then present analytic estimates and numerical 
calculations of the returning positron fraction and PC heating luminosity for CR screening and 
compare our predictions with ROSAT observations. 

The most recent relevant study is that of Hibschman & Arons (2001), who focus on the 
determination of theoretical radio pulsar death lines. Even though they use the frame-dragging 
electric field in their treatment of primary electron acceleration and incorporate both CR and ICS 
photons in the pair production, our approaches are essentially different. For example, in their 
calculation of the PC heating they use a very rough estimate of the returning positron fraction 
(which is generally inaccurate and is the same as that given by MT92). Also, they refer to the 
PFF as the location where the electric field is fully screened, whereas in this paper we define 
the PFF as the location where pair production begins. In the case of CR, this distinction is 
minor since the screening scale is small compared to the PFF height, but in the case of ICS the 
distinction is very important since the screening scale is comparable to the PFF height. Thus, the 
quantitative results of our studies are not directly comparable, although some of their conclusions 
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are qualitatively similar. A more detailed comparison of our results with theirs will be given in 
our next paper, where we will present our theoretical radio pulsar deathlines. 

In Section 2, we describe the model we employ to calculate the accelerating electric field in the 
PC region of pulsars. We outline our main assumptions and approximations. The electrodynamic 
boundary conditions at the stellar surface and upper PFF are specified. In Section 3, we present 
our calculations of the returning positron fraction. In this paper our analysis is limited to 
the parameter space where the CR is the primary source of photons producing pairs in the 
magnetic field. We give our analytic estimates first. Then, we discuss our numerical treatment 
of electron-positron pair dynamics and effect of partial backfiow of positrons. In Section 4, we 
calculate both analytically and numerically the thermal X-ray luminosities from pulsar PCs 
produced by precipitating positrons. Summary and conclusions are given in Section 5. 



2. THE UNDERLYING ACCELERATION MODEL 

2.1. Electric fields in tiie regime of space-charge limitation 

In this Section we shall summarize some of the general solutions pertaining to the 
electrodynamics of open field line regions of rotating NSs derived earlier (MT92, MII97, and 
HM98) and which will be exploited throughout this paper. We should note that these solutions 
are based on the assumption that a regime known as space- charge limitation of current occurs in 
the acceleration of electrons ejected from the PCs of a NS. The effect of space-charge limitation 
is well-known since the very first experimental and theoretical studies of the vacuum diodes, and 
detailed discussion of this effect can be found in any textbook on electronic devices. In Section 2.2 
we will briefiy outline some of the features of the space-charge limitation specific to the physical 
conditions of the pulsar PC regions. 

The general relativistic Poisson's equation describing the electric potential distribution (in 
the corotating frame of reference) within the region of open field lines in the magnetosphere of a 
rotating NS reads (MT92) 

V- (^ivcD^ =-47r(p-pGj), (1) 

where p is the space charge density, pcj is the general relativistic analog of the Goldreich-Julian 
(GJ) space charge density, a = y^l — e/ry is the redshift function, e = rg/R, rg is the gravitational 
radius of the NS, R is the stellar radius, rj = r/R is the dimensionless radial coordinate, 
K = el /MR?, I and M are the moment of inertia and mass of the NS, respectively. The differential 
operators such as gradient, divergence, etc. should be taken in corresponding curvilinear 
coordinates. Eqn. (|l|) takes into account both the effect of dragging of inertial frames of reference 
and gravitational redshift. The former dramatically affects the electrodynamics of a NS, while the 
latter rather moderately modifies the corresponding electrodynamic quantities. 

The study of particle acceleration within the region of open magnetic field lines in pulsars 
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is grossly facilitated by the fact that for all known pulsars (including ms pulsars) the angular 
size of the PC is less than ~ 0.2/-\/ P/3 ms radian (P is a pulsar spin period). Thus, we can 
use a so-called small-angle approximation, which proves to be very satisfactory for most relevant 
problems. In this approximation, for general relativistic dipolar magnetic field, 
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where Vt is the angular velocity of the NS rotation, Bq is the surface value of the magnetic field 
strength at the magnetic pole and x is the pulsar inclination angle. Here we use the dimensionless 
radial coordinate r] (= r/R) and characteristics describing field-streamlines ^ = 9/9{rj) [which 
is a magnetic colatitude scaled by the half-opening angle of the polar magnetic flux tube, 
9{7]) = Oo^/r|f{l)/f{r]), where = \/^B,/cf{l) is PC half-angle], (p is magnetic azimuthal angle, 
and the functions f{rj) and H(r]) are factors accounting for the static part of the curved space-time 
metric [see Eqns (All), (A12) in Appendix]. 

In this paper we continue our study of a regime of space- charge limitation of current, which 
allows us to derive an explicit expression for a steady-state distribution of space charge density 
(see MH97 for details) 
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Using expressions (^)-(^ and appropriate boundary conditions we can solve Eqn. (|^). In this 
paper we will be mostly interested in the solutions for the accelerating electrostatic potential 
(field) valid for altitudes less than one stellar radius from the PC surface i.e. for z = rj — 1 < 1, 
even though their generalization for z > 1 proves trivial (HM98). Thus, for z < 1 the general 
solution for the potential can be expressed as 
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Here <I>o = {^R/c)BqR is the dimensional factor which is nothing but an order-of-magnitude 
maximum value of electrostatic potential generated by a magnetized globe (with dipolar external 
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magnetic field) rotating in vacuo (see Deutsch 1955); ki and ki are the positive zeros of the Bessel 
functions Jq and Ji, respectively, indexed in ascending order; 6{r]) is another correction factor [see 
Eqn. (A13)] accounting for the gravitational redshift effect, with ^^(1)15(1) 1 (see e.g. HM98); 
and Oq is defined right after Eqn. (^). The explicit expressions for the function J^i are determined 
by the solution of Eqn. (|l|) and appropriate electrodynamic boundary conditions. The solutions 
for the function J^i relevant to our study will be presented in § 2.3 and the Appendix. In this 
paper we do not need and therefore do not present the explicit general-relativistic expressions for 
the magnetic field, but we shall refer an interested reader to our previous publications (MH97, 
HM98). Finally, our results will be based on the electrodynamic solutions of Eqn. ([l|) subject to 
the standard Dirichlet boundary conditions. 



2.2. Modification of space-charge limitation by returning positrons 

In this Section we shall discuss the space-charge limitation of current taking into account 
the returning positron flux from the upper PFF. The idea of space-charge limitation of flow 
was in fact introduced into pulsar physics by Sturrock (1971): " What happens to the positrons 
produced by the electron-positron avalanche? If they were all returned to the surface, as one might 
expect by consideration of the electric field, the resulting space charge would reverse the sign of 
the electric field at the surface, cutting off the flux of primary electrons". Later on this idea was 
investigated by Michel (1974, 1975), Tademaru (1973), and Cheng and Ruderman (1977), and 
then substantially quantifled in the calculation of the electric fields produced by the relativistic 
electron beam within the region of open field lines in a NS magnetosphere (Arons & Scharlemann, 
1978; Arons 1983, hereafter A83). They found that only a small fraction of positrons (e.g. ~ 6q 
of the negative charge density already present in the electron beam, as estimated by Arons & 
Scharlemann) are returned relative to the primary electron beam. Note that the calculation of the 
fraction of positrons returning from the upper PFF is far from being a trivial problem. In this 
paper we attempt to address this problem again, taking into account most recent updates of the 
physics involved and revising some of the underlying assumptions. 

As a result of injection of backflowing positrons into the initial electron beam, an excessive 
negative space charge unavoidably builds up which lowers the electrostatic potential above the PC 
surface. The lowering of potential above the surface and forcing its gradient to reverse sign at the 
very surface would inhibit electron emission, thus resetting the electric held to zero at the surface. 
The fundamental difference between the laboratory vacuum diode and PC voltage generator is 
that in the former the thermionic cathode is the only emitter of electrons, whereas in the latter 
the PC surface, serving as a cathode, both emits electrons and collects positrons from the upper 
PFF. Thus, for a vacuum diode the space-charge limitation of current is solely determined by the 
cathode, while for a NS PC this regime of operation depends on the upper PFF as well. Because 
of the intrinsic limitation of the total current by the GJ value, the injection of counterstreaming 
positrons suppresses the primary electron emission: the primary electrons should be ejected from 
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the stellar surface with a density somewhat less than the GJ charge density to keep -Ey = at the 
surface. This effect provides a negative feedback between the primary electron ejection and flux 
of returning positrons, since the primary electrons ultimately determine (through the electrostatic 
potential drop they generate above the PC surface) the fraction of returning positrons. 

Now, using very basic arguments, we would like to demonstrate how the accelerating electric 
field is related with the fraction of returning positrons. In our previous paper (see HMOS) we 
quantitatively explored the effect of rescaling of the lower boundary (e.g. stellar surface) caused 
by the formation of a lower pair front in the regime of space-charge limitation. The returning 
positrons add to the space charge density of primary electrons and lift up the boundary at which 
the p = pq:] condition is satisfied, similar to the effect produced by the formation of lower pair 
front. Let us use this reasoning and express the space charge density of a primary beam in terms 
of effective (rescaled) stellar radius, Re- For simplicity, we illustrate this by employing solution 
(^ at X ~ 0- The expression for the space charge density at the rescaled stellar surface then 
reduces to 

p = -a{l-r,lK), (9) 

where rj^, = R/Re (< !)• 

If there were no returning positrons (77* = 1) we would simply have 

p_ = -cj(l - k). (10) 

The returning positrons perturb the primary electron beam causing extraction of less electrons 
than in the case of unperturbed electron beam and effectively mimic electrons missing in a beam. 
Assuming that the regime of space-charge limitation is maintained, the following condition should 
be satisfied at the surface 

p = p_+5p+, (11) 

where 5p+ is a small perturbation of (and of the same sign as) the electron space charge density 
caused by returning positrons. Introducing the dimensionless fraction of returning positrons, x+, 
so that (5/9+ = — (7KX_|_, from the above equation we get that rj^, = (1 — x+)^/^. Thus the positrons 
returning from the PFF reduce the accelerating electrostatic potential (field), so that the latter 
gets factorized by (see also HM98) 

= 1 - (12) 

It is now clear that the net effect of the injection (e.g. from the upper PFF) of counterstreaming 
positrons into primary electron beam is a weakening of the source of effective charge, or, 
equivalently, the transforming of general-relativistic parameter k into (1 — x+)k. Note that x+ <C 1 
in most cases. Thus the final result does not depend on the particular value of rescaled stellar 
radius and illustrates the dynamics of the feedback between accelerating electric field and fraction 
of returning positrons. Namely, the returning positrons tend to reduce the accelerating electric 
field and therefore suppress the ejection of primary electrons. The starving electron ejection is 
compensated by returning positrons, which mimic the outflowing electrons, thus recovering the 
original potential drop. Obviously, the increase in the fraction of backfiowing positrons reduces 
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the efficiency of pair production by primary electrons wfiich may eventually cut off tfie supply 

of positrons themselves. One can expect some demand-supply balance to be established in the 
beam between the densities of primary electrons and backflowing positrons. Thus, we are now 
in a position to calculate acceleration of primary electrons and returning positrons produced 
at the PFF in a self-consistent way, taking into account the weakening of the primary beam 
by returning positrons. It is important that, after the appropriate parameters that favour the 
steady-state regime are found, we can use essentially the same expressions for the accelerating 
electric field/potential as in the case of a pure electron beam, simply because in a steady-state 
electrodynamic solution the returning positrons arc undistinguishablc from outflowing electrons. 
Then, the problem reduces to finding (through numerical iterations) a stable solution (see § 3.2.2 
for details). 

2.3. Screening condition at the upper PFF and accelerating electric field 

In this paper we derive an appropriate solution for the accelerating electric field, and use 
this solution for a self-consistent analysis of a steady-state flow of charged particles (electrons 
and positrons) from the PC surface up through the PFF. We assume that both the PC surface 
and surface formed by the last open magnetic field lines are equipotcntial ($ = 0). We must 
note that, within the context of the lower boundary condition, the PC surface is actually either 
the upper boundary of the stellar atmosphere or the bare stellar crust. Since the free emission 
of charges from the PC seems to be more favourable than charge starvation above the PC, the 
electron current should be consistent with the global magnetospheric current. We assume that in a 
steady-state regime the rate of ejection of electrons from the PC surface is approximately equal to 
the rate of inflow of electrons from beneath the surface. It is this basic condition that makes the 
ejection of electrons occur without producing significant polarization of charges at the PC surface. 
Note that, in a static model ignoring the global current closure the free emission of electrons from 
the PC would rather result in a filling of the polar magnetic flux tube with electron-ion plasma 
and any electron acceleration will be halted. 

The upper boundary condition also needs to be discussed in detail, simply because it proves 
to be crucial for the calculation of the fraction of returning positrons from the upper PFFs. This 
condition was extensively discussed in the literature (see e.g. Arons & Scharlemann, 1979; A83) 
and it was suggested to set = and Vy • i^y = at the upper PFF. The problem is that these 
two conditions arc formally inconsistent, as we shall demonstrate in the next paragraph. This 
inconsistency should not be ignored in any more or less realistic model simply because these two 
conditions cannot be satisfied not only at the same point but even at two points separated by a 
distance of order of the PFF thickness. 

The formal condition assumed previously that p = poj at the upper PFF (at r] = rjo) means 
that in the very vicinity of the PFF the radial part of electrostatic potential should satisfy the 
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following equation [see e.g. MH97, Eqn. (48)] 



d'^P dP n 

—^ + 2—-a^P = 0, 13 
a?]^ or] 



where 2/a = 0(?7o)%-\/l — e/% ~ I'pc/ R i^pc is the PC radius). The non-trivial solution of this 
equation reads 

P = Cexp{-r]) sinh[VTTa^(??o - ??)], (14) 

where C = const. The boundary conditions = and Vy • -By = at ?] = r/o are simply reduced 
to the following conditions for the function P 

dP , d'^P dP , , 

— = and — - + 2— = 0, 15 
df] df]'^ drj 

respectively. One can easily verify that the first of the above conditions is not satisfied. This 
means that it is incorrect to assume that both Vy ■ E^^ =0 and -By = are satisfied at the same 
radial coordinate. The reason is that condition Vy • =0 actually implies by Eqn. ( p!^ ) that the 
electrostatic potential also vanishes at the upper PFF. However, the latter by no means requires 
that vanishes at the same boundary. Moreover, as it is illustrated in Appendix [see Eqns (A2), 
(A3), and (A7)-(A10)], the vanishing of both Vy • i?y and i?y is very unlikely even within entire 
PFF. In general, specifying both Dirichlet and Newmann conditions (e.g. Vy • £^y =0 and -Ey = 0) 
overdetermines the problem and leads to there being no solution. 

In our solution (which is essential for the kind of problem we discuss in this paper) we will 
only require that beyond the upper PFF, 

<^ ^ <I>oo = const / 0, as P^PgJ) (16) 

where $00 is the potential at altitudes well above the upper PFF. This condition implies that 
there is no perfect adjustment of the effective space charge density to the GJ space charge density 
at the very onset of the upper PFF. When i?y has dropped to a low enough value that positrons 
are no longer able to turn around, then there is still non-zero ^ but no additional returning 
flux. If, at the PFF, the GJ space charge were exactly compensated by the effective space charge 
of primary electrons and electron-positron plasma, then the electrostatic potential would drop 
to zero at the PFF. The vanishing of potential at the PFF suppresses the accelerating electric 
field and even causes it to reverse its sign well below the PFF. The grounding of potential at 
the PFF would therefore dramatically affect the acceleration of primary electrons: the electron 
beam will be cut off from the PFF by a layer with a reverse polarity of the electric field, and 
the primary electrons will no longer be able to maintain the PFF they would produce if there 
was no strong negative feedback under discussion. As a result, it is very unlikely that a stable 
upper PFF would be established at all, if the condition p = pQj sets in at the PFF. We suggest 
that requirement (^) is more favorable for a non-disruptive regime of acceleration of primary 
electrons and formation of upper pair fronts. This requirement allows the parallel component 
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of the electric field to partially penetrate into the relativistically moving pair plasma, with the 
bulk of the electron-positron pairs streaming out as a quasi-neutral beam. We therefore relax the 
intuitive requirement that relativistically moving electron-positron pairs completely screen out the 
electric field at or right above the PFF. A self-consistent treatment of the dynamics of primary 
electrons and pairs in a nearly screened electric field, as will be discussed below, may provide us 
with a reasonable estimate of the fraction of returning positrons. 

Using general expression ^ we can derive the formulae for the accelerating potential (field) 
for < z < 1. This solution has the property that -Ey saturates (-By const) for rpc/R « z < 1, 
before it declines approximately as (see MT92, MH97, HM98). In this case (no explicit upper 
boundary) the solution for the function J^i(z,ji) is 



J^i = liz + exp(-7j2;) - 1, 
and the corresponding formulae for the potential and field read 
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where Eq = VlRBq/c = ^q/R. 



(19) 



Simple expressions can be derived from the above formulae in some limiting cases. For 
example, for 72; ^ 1 (z <C rpc/R) we get 
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In derivation of the above expressions we have used the formulae 
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fitting the results of numerical summation with accuracy better than 2%. For z ^ Tpc/R (72; ^ 1) 
we get 



3 Oi? z 
2 °~7(1) 
3 J^iZ 1 



fi 



(i-a 



K cos X + ^^o^^(l)<^(l)^ sin X cos ( 
K cos X + 76'o-H'(l)'^(l)^ sin x cos ( 



(24) 



(25) 



Before we proceed with our modeling of field screening and calculation of the fraction of returning 
positrons, let us introduce a dimensionless minimum height zq (= 770 — 1) at which a first pair is 
produced . Then for our modeling we will need the expression for the electric field applicable to 
the region of pair formation. Since the penetration depth of the electric field into the pair plasma 
(or the characteristic length scale of the electric field screening) is much less than the characteristic 
length scale of the accelerating electric field [e.g. Eqn. (]25|)], the electric field in the pair region 
can be satisfactorily described by the formula 



Ef{z ^ zo) = (zo) exp[-(z - zo)/A.], 



where E^^^^{zq) is the accelerating electric field given by ( [l^ ) and calculated at z 
characteristic length scale of field screening. 



(26) 

zq, As is the 



3. CALCULATION OF RETURNING POSITRON FRACTION 

3.1. Analytic estimate 

Let us estimate very roughly the fractional density of positrons flowing back from the PFF. 
We shall restrict ourselves to the case where the pair creation is mostly determined by the CR 
photons. The condition that the positrons with energy e+ (in mc^) turn around within the PFF 
can be written as 

e\Eii\As = e+mc^ , (27) 

where i^y is the accelerating electric field evaluated at the PFF, is the same as defined in Eqn. 
(p^, and e = |e| is the elementary electron charge. The characteristic energy of turning around 
positrons can be estimated as 

2 n+ aNe 

where is the characteristic energy of pair-producing curvature photons, rig is the number density 
of primary electrons, n_|_ is the number density of turning around positrons, and dN^/dNf. is the 
number of photons (per primary electron) producing those pairs whose positrons turn around and 
flow back to the surface. Note that 

dN^ , , Ae^ , , , 
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where p{s) = {V^e^ /2Trpc)jF{£/ecr) is the spectral power of CR (see e.g. Landau & Lifshitz, 1975), 
Ae^ is the characteristic interval of energies in the CR spectrum from which the pair-producing 
photons generate positrons with the characteristic energy ~ e+. Using Eqns (p7|)-(p9[) we can write 



which translates into 



n+ 1 Ae-^ mc 



^ « 5 X W~'AF{e,/e,r) ) , (31) 



where Scr = •i^cl^/'^Pc (^c = h/mc = 3.9 x 10 cm is the Compton wavelength) is the critical 
energy of the curvature spectrum. The factor A (= Ae-y/scr) in Eqn. (31) takes into account 



the fact that only photons with approximately Gaussian distribution around some characteristic 
energy effectively produce pairs, with Ae^ being the spectral interval of photons that produce 
returning positrons. The factor F{x) [= x K^/^{z)dz vvrx/S exp(— x), if x » 1, where 
K^/^ is the modified Bessel function of order 5/3] accounts of the fact that returning positrons are 
produced by CR photons with ciiGrgies greater than Scr- 

Finally, factor 2e7'^/3/9^|£'|| | < 1 describes 
the efficiency of emission of CR by accelerating electron. It reaches a maximum in the so-called 
"saturation" regime where the power of electrostatic acceleration is almost exactly compensated 
by the CR losses. If we adopt in Eqn. (|3l| ) A 3, and ~ (3 — 6)ecr (this range for of photons 
producing returning positrons agrees with that resulting from our numerical simulations, with the 
characteristic value of e-y being a factor of 2-3 smaller than the mean energy of pair-producing 
photons given by Eqn. [29] in HM98), we get 

^ « (0.001 - 0.02) (,^^] . (32) 

It is interesting, that the fractional density of returning positrons given by the above rough formula 
is practically independent of pulsar parameters, except for the dimensionless factor in brackets. 
The latter also tends to be close to its maximum value of unity. The reason is that in most 
physically interesting situations 7 should be well tuned to allow magnetic pair production, which 
occurs at or near the saturation part of the acceleration curve, where the further acceleration of 
electrons would be suppressed by the CR losses, i.e. where the condition e\E^^\ ~ 2e^7'^/3p^ is 
roughly satisfied. 

Even though formula (^) gives the right range for the fractional density of returning 
positrons, it is based on a simplified mapping between the spectrum of pair-producing CR photons 
and the distribution function of returning positrons. Its main deficiency is that it requires a priori 
knowledge (say, based on numerical simulation) of parameters A and e^. Also, it does not reflect 
the fact that at the upper PFF the maximum density of returning positrons is limited by the 
value of \pGj — p\ and is therefore dependent on the altitude of the PFF. In a steady state one 
can estimate the maximum value of the density of returning positrons as one half of the difference 
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\pgj{zo) — p{zo)\- Thus, using expressions (|2|) and @ for pcj and p we can come up with an 
alternative formula 



Pgj(2o) 2 



2 1 



—zo. (33) 



A similar formula has been used by Zhang &; Harding (2000) in their study of the X-ray 
luminosities from the spinning-down pulsars. 



Let us perform numerical estimates of zq based on formulae (21) and (^5|) for the accelerating 
electric field. To make the resulting estimates as compact as possible, we shall adopt the following 
parameters: ^ = 0.5, cosx = Ij and k = 0.15. This is still a justified simplification, since the 
resultant expressions will have rather weak dependence on these parameters. Then, Eqns (^ij) and 
( p5D reduce to 

i?l|6 = 42^z, (34) 

-^0.1 

and ^ 

% = 0-57#' (35) 

respectively Here E^q = E^i/W^ esu, Bu = Bq/W^"^ G, and Pq.i = ^/O.l s. 

The altitude of the PFF above the stellar surface can be estimated as [see e.g. HM98, Eqn. 

(1)] 

So = min[Sa{jmin) + Sp{emin)], (36) 

where Saijmin) is the distance required to accelerate the particle until it can emit a photon of 
energy Smin and Sp{emin) is the photon pair-attenuation length. Thus, using expressions (|3^, 
( pBD we arrive at 

So = min{Sa + B^/j'^), (37) 

where 

ci,ii = / ^7'/' if Ell is given by (21), 

\ A}l-f if is given by (H). ^ ' 

Here A\ = 9.1 • Pq^.i^/sJ^ , Al^ = 3.3 • W^^PlJBu, and B^ = 1.95 • W'^^Po.i/Bu. Note that this 
expression for B^ is valid for B12 < 4.4. Equation (^) has a minimum at the value of Lorentz 
factor 

\ 2.1 Pj-;'" it E|| is given by (||). ^ ' 

By substituting expressions for Sa and 7mm into (^) we get 



1.9 Poy^^^i//^ if E\\ is given by (21 
3.0 Po{^B^2^ if is given by (H 



.0 = 5o/i2 = 0.03 ^ - -ji :^ (4o) 



Note that formulae ( ^I[) and (25) for the accelerating electric field are most likely applicable when 
zq < Tpc/R and zq > rpc/R, respectively. Using the above expressions for zq, these two criteria 
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can be reduced to -Pq < 0.5 Bu and Pq^i > 0.4 B12, correspondingly. Then, using relation 
P0.1/B12 ~ 1.2A^/tq {tq = r/lO^ yr, r = P/2P is the pulsar spin-down age) we can rewrite formula 
® as 

'"^"•"'"^ il.lPo'.f if Po^.f > 0.4 i^,,. 

For a Crab-like PSR with r ~ 10'^ yr, Pq.i = 0.3 and P12 = 8 the above formula yields 

zq k 0.007, (42) 

whereas for middle-aged PSRs with r ~ 5 • 10^ — 10^ yr, Pq.i =2 — 3, and B12 = 7 — 8 we get 

zo ~ 0.01 - 0.02. (43) 

Finally, for an old 10-ms PSR with r ~ 10^ yr and B12 = 8 • 10"^, we get 

ZQ ^ 0.2. (44) 



Inserting Eqn. (41) into Eqn. (^) we arrive at the following expression for fractional 



returning positron density 



~ 01rV2 / 1-85 bU'P.T' if Pl{' < 0.5 Bu, . . 

^r'-''^' \3.0P^{' if Pli' > OA B,2. ^ ^ 



3.2. Numerical calculation of the returning positron fraction 

3.2.1. Pair source function in screening region 

The pair source function Q^{'yQ' ,Xnr) is the joint initial energy and spatial distribution 
of electron-positron pairs produced by the primary electrons, where Xnr is the distance above 
Zq R. Only the first generation of pairs is important to the screening process, as the higher 
pair generations are produced beyond a distance of zq + A^^^A^ [Ng is the number of screening 
scales specified in the end of Section 3.2.2). It can be shown that the attenuation length of 
synchrotron photons from the first pair generation is much larger than A^. Each primary electron 
is accelerated from its starting point at the stellar surface, radiates CR photons and the pair 
production attenuation lengths L of these photons are computed, as described in detail by HM98. 
The location of the first pair defines the PFF at zq. We compute the pair source function by 
accumulating the number of pairs per primary particle at each energy and height above the PFF 
as the particles accelerate up to the PFF, averaged over simulations for 10 - 50 primary electrons. 
The primary particle trajectory is divided into discrete steps and at each step the CR spectrum 
at energy 7(5) is divided into discrete energy bins. A representative photon from each CR energy 
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bin is propagated through the local field to determine whether it produces a pair or escapes (for 
details of such a calculation, see Harding et al. 1997). We accumulate a survival probability, 



where 



Psmv{s) = exp|-T^(s) j, (46) 



r(0kB,c^)ds' (47) 



is the optical depth along the path. The PC angles of ms pulsars are large enough that the 
curvature of the photon trajectories in the strong gravitational field of the NS is important and 
we take this into account in computing the pair attentuation lengths (see Harding et al. 1997, 
for details). A random number ?R. is chosen to determine the pair production point of each test 
photon, when Psurv(s = L) = ^. The momentum of each pair member is assumed to have half the 
energy and the same momentum, parallel to the local magnetic field, as the parent photon. The 
"number" of CR photons, n-cR, represented by the test photon in each energy bin is estimated by 
dividing the energy radiated in the spectral interval, Ae* = eJnax ~ ^min' by the average energy in 
that interval, (e*), 

-Y As fi^^'' ^CR(.e)de 
(e*)c Jo NcR{e)de 

where 

NcR{e) = - jFie/scr) (49) 

c 

is the CR energy spectrum, Scr is the critical energy and pc the field line radius of curvature 
[see also the formulae right after Eqn. (31) ] and is the CR loss rate ^f^^mc^ = 2e^c7^/3/3^. 
The height above the PFF and the energy of the pairs are accumulated in a two-dimensional 
distribution, normalized by the total number of test primary photons and by ncR, to form a 
distribution Q^{'y^,Xnr) of the number of pairs per energy per primary at each height and energy 
interval. An example of one of the computed pair source functions is shown in Figure 1. The first 
pairs have the highest energies because they have the shortest attenuation lengths. The number 
of pairs increases rapidly beyond the PFF and the mean energy of the pairs decreases. 



3.2.2. Pair dynamics 

The charge density along open field lines above the PFF will increase due to injection of pairs. 
The dynamical response of the pairs to the -Bjj'^, i.e. acceleration of electrons and deceleration of 
positrons, will determine the increase in effective space-charge density with height in the screening 
region. We find that the E^i^ is large enough that the major contribution to the space-charge 
density results from positrons that are turned around and accelerated back toward the NS surface. 
Above the PFF, the continuity equations for elections and positrons are 
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dn dn d{n (3 c) 

where n'^ and n~ are the first moments and and /5~ are the second moments of the positron 
and electron distribution functions, 7^(7^, x), respectively: 

P^ = — =[ -f^{l^,x)dv (52) 
c J c 

n^ = I f^i7^,x)dv. (53) 

Since the source function is electron positron pairs, n"*" = n~, we can subtract Eqn. (^) from 
Eqn. ( [50| ) to give 

^d{n+ -n-) _dp dn+P+ 9n-/?-\ 

dx dx \ dx dx I ^ 

so that the total charge density is 

p = e{n+P+ -n-p-). (55) 

The distribution functions of electrons and positrons f^{v^,x) are computed by dynamically 
evolving the pair source function, subject only to the force of the electric field in the screening 
region. Pairs are injected at discrete points, Xnr, and with a discrete distribution of energies, 7^ 
as described below. Each sign of charge in each of the Ug energy bins is evolved in energy through 
a separate grid, x = {z — zq)R, of E^^{x) computed at discrete points x„. Using Eqn. ([2^), we find 
that with the resolution of the calculation, the region in which the positrons are non-relativistic 
at their turn-around points is not resolved. We are therefore justified in treating the particles as 
relativistic everywhere and may thus ignore the momentum equation. The energy of a particle at 
point x in the grid is 

7^(x) =lo^^ r El'^{x')dx', (56) 

so that 



1/2 

1 - 77^^ \ ■ (57) 



V (x) 

Particles from the source function, Q^{'y^ evolved through the grid until they either 

travel upward across the top boundary or downward across the lower boundary. At each grid 
height, the steady-state distribution function of electrons and positrons are computed 



/-(7-,x) = ^^^||l^. (58) 



nr 



Using Eqns (p^), (|53|), (|5q ) and (|58|), we can determine the density of upward-moving electrons, 
n~(x), and positrons, n^{x) and the downward-moving positrons, f^^(x) at each grid point. The 
total charge density due to pairs is then 

p{x) = e[n+(x)/3+(x) - n~{x)P~{x) + nj{x)(3^{x)]. (59) 
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The downward-moving positrons will constitute an additional upward-moving negative 
current which will add to the current of primary elections (Section 2.2). Because the net current j 
is required to be the GJ current, we therefore need to adjust the primary current for the returning 
positron current. Thus, we can write 

n;(l)=n;-n+(l), (60) 

where n~ is the number density of primary electrons and ra^(l) is the downward- moving positron 
density at the grid lower boundary (the PFF). The charge deficit, Ap = Pp — pcj at the PFF is 
thus also readjusted. 

The calculation thus proceeds as follows. We first compute the location of the PFF above 



the surface, zq, and the pair source functions as described in Section 3.2.1 . An initial screening 
scale length, A^, is chosen as a fraction of zq, the distance to the PFF. The electric field in the 
screening region, above the PFF, is modeled by an exponential (Eqn. |^6|) with scale length A^. 
The pair dynamics determines a pair distribution function and a returning positron density, as 
described above. The primary electron flux and Ap are adjusted for this returning positron flux, 
and the pair distribution function and a returning positron density are then recomputed using 
this adjusted charge deficit. A new screening scale length is determined from the point where 
p{NsAs) = Ap, where Ng is the number of -Bj|'^(x) scale heights required to guarantee that no 
positrons turn around above zq + NsAs- The solution will then be self-consistent. We find that 
Ns = 4.0 satisfies this condition and set it as a constant value. The scale height of the electric 
field in the screening region is then set to the new A^. The iteration continues until the returning 
positron density and screening scale attain stable values. 



3.2.3. Screening scale height and returning positron density 

Convergence of the screening scale height A^, the total charge density above the PFF, 
p{x)/pG.j, and the relative returning positron density, p^{l), is achieved independently at each 
colatitude, ^, across the PC. An example of a self-consistent solution for p{x)/ pcj and A^ is 
shown in Figure 2. The charge density increases rapidly immediately above the PFF, where the 
pair source function is growing exponentially and where E^^^ is high enough to decelerate and turn 
around all the produced positrons. At a distance above the PFF comparable to the screening 
scale height, E^^ has decreased to a fraction of its value at the PFF and the rate of increase 
of the charge density begins to moderate, as fewer positrons are able to turn around. We find 
that the creation of the charge density and subsequent screening of E^^ is due almost entirely to 
returned positrons rather than to velocity differences between electrons and positrons as the pairs 
dynamically respond to the electric field above the PFF. Due to the boundary condition imposed 
on the potential (see Eqn. [|l6|), E^^ never decreases to 0, but only approaches at infinity. 
Consequently, the self-consistent charge density never achieves the value Ap but approaches it 
asympotically from below. 
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Figure 3 shows solutions for the screening scale, A^, and the returning positron density, 
p+/ PGJ-, across the PC as a function of the colatitude ^, scaled to the PC opening angle at the 
surface, for various values of the pulsar period and surface magnetic field strength. Both and 
p+/ Pgj increase toward ^ = (magnetic pole) and ^ = 1 (the outer edge of the PC). Since the 
value of the GJ density, pcj is constant with the variation represents a true variation in />+ 
across the PC. The increase in A^ and p+ toward ^ = is due to the increasing field line radius 
of curvature near the pole, causing the pair attenuation length to grow large. The pair source 
function grows more slowly with distance above the PFF, and the screening scale height increases. 
The returning positron density increases because as the height of the PFF increases, the charge 
deficit also increases, requiring more returning positrons for screening. Near the pole, there are no 
solutions for A^ and /)+ and screening is either incomplete or there is no screening at all. The 
increase in Ag and /9+ toward = 1 is due to the decrease in E\\ caused by the boundary condition 
$ = at ^ = 1 imposed on the solution to Poisson's equation. Consequently, the primary particles 
must accelerate over a longer distance to produce pairs, increasing the height of the PFF and 
thus /?+. Also the positron turn-around distance becomes longer with a lower E\\, increasing the 
screening scale height. 

However, the variation in A^ and p+{C) with ^ is small compared to their variation with 
pulsar period and surface field strength. Comparing Figures 3b, 3c, and 3d, one can see that the 
screening scale height A^ increases by a factor of about ten as the period increases by only a 
factor of four, for a constant surface field. Comparing Figures 3a and 3c at a period of 0.1 s, Ag 
decreases by a factor of about six as the surface field increases by only a factor of two. Figure 4 
shows the dependence of the screening scale length, A^, on surface field and period. Generally, 
As decreases with increasing surface field strength because the pair attenuation length is shorter 
in higher fields. As a result, the pair density grows faster and E\\ can be screened in a shorter 
distance above the PFF. Consequently, ms pulsars have relatively large screening scale lengths, 
but As is still small compared to zq. Figure 5 shows the dependence of p+ on surface field and 
period. While is less sensitive to pulsar parameters than A^, it decreases significantly as the 
surface field increases and as the period decreases. The numerical values of p+Z/OGJ are smaller 
than, but within a factor of two of, the analytic estimate of Eqn. (^). The variation of A^ and 
p+ with pulsar parameters is primarily due to the dependence of E\\ and the pair production 
attenuation length on period and surface field strength. By contrast, Ag and p+{Cj vary by only 
factors of two or three across the PC of a single pulsar. 

Since />+ is always a small fraction of both the primary particle density and of pcj) only the 
very first pairs of the full, multi-generation cascade are needed to screen the E'n . We therefore did 
not simulate the full cascade in computing the pair source functions needed in the present work. 
The vast majority of the cascade pairs are produced in the region where E\\ ~ and they freely 
escape the magnetosphere with only radiation losses, but no significant acceleration. Previous 
simulations of the full cascade produced by the primary particles (Daugherty & Harding 1996) give 
multiplicities ~ 10^ — 10^ pairs/primary. The fraction of the total number of secondary pairs that 
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are returned to heat the stellar surface can thus be estimated by dividing the values of /9+/pGJ in 
Figure 5 by the cascade multiplicity. 



4. POLAR CAP HEATING RATES AND THERMAL X-RAY LUMINOSITY 

4.1. Analytic estimate 

Let us derive an expression for the total power that can be deposited onto a single PC by 
precipitating positrons. The general expression can be written as 



ac [ p+izo,OHzo,OdS, (61) 



where the integration is over the area of a sphere cut by the polar flux tube at the radial distance 
rjQ, and the factor a accounts for the general-relativistic correction to the current (j ~ acp). Here 
dS = [Spcr]Qf{l)/7rf{r]o)]dQ^, Spc = ir^R^ /cf{l) is the area of the PC, dQ^ = ^d^dcj) is an element 
of a solid angle in the PC region. Thus expression (|6l] ) reduces to 



L+ = 2acSpcVo7r\ t ^+(^0,0^(^0, Oidt (62) 
/(%) -'0 



After inserting expression for $ [see Eqns (|20|), (|24|)] into (^) and normalizing by pcj, we get 

L+ = f+Erot, (63) 



where 



„ /»,3^^2 P+{zg,C) A2.19^0.705( 



and 



U = —j^^{l-^l4)4^os^X ^-^^{l-e-^y-'-'idi (64) 
^oVl-e "'0 PgAzo) 

U = 9^(1 - ^^HDzo cos^ X t P±^^{1 - em (65) 

are the fractions of pulsar spin-down power consumed by returning positrons, corresponding to 
the cases where $ is given by Eqn. ( |20| ) and (|2^), respectively. Here 

is the general-relativistic expression for the pulsar spin-down losses in vacuum (see MH97). In 
this formula Bq/ f{\) is the surface value of the magnetic field strength "as seen" by an infinitely 
remote observer, whereas Bq is the value measured locally, at the NS surface. In a flat-space limit 
Bq/ f{\) simply transforms into Bq, and we get a classical formula for the magneto-dipole losses. 
Note that in our derivation of formulae ( |64|) and (^) we used only the components of potential 
and charge density that are proportional to cos x and we assumed that zq is independent of ^ (flat 
PFFs), which is still a satisfactory approximation for rough analytic estimates. Our numerical 
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calculations of /+ take into account the curvature of PFFs (see further discussion of the numerical 
results), though. 



Now, using expression (^T|) for zq and ( |45| ) for p^/pcj) we can combine Eqns ( |6^ ) and (65) 
and write (setting cos x — 1, — 0.15) 

^ 3 / 0.88 (Po^/re)^/^ if ^af < 0.5 P12, 

This expression illustrates that there is a rather strong correlation between the efficiencies of PC 
heating by precipitating positrons and pulsar spin-down age. The important implication of this 
correlation is that older pulsars should favour enhanced PC heating (see, however, the warning 
on use of Eqn. |6^ in Section 4.2). It is interesting that the above order-of-magnitude analytic 



estimates are in good agreement with our numerical calculation. 

A81 derived an expression for PC heating luminosity, taking into account the curvature of 
field lines but not general-relativistic effects in the derivation of -By. Using his results we can 
derive the following expressions for /-|_ 

The value of in Arons' model also increases with pulsar age but is significantly lower than our 
result. This is because the frame dragging electric field and consequent acceleration energy in our 
model is higher. 

During the photon cooling era, the effective temperature of a NS decreases with age according 
to a power law. For example, using the cooling calculations by Page and Sarmiento (1996), we can 
estimate that the cooling luminosity, Lcooi, roughly scales as Tg~^ (rg > 0.3) and Tq ^ (rg > 1) for 
the case with and without core neutron '^P2 pairing, respectively. By normalizing the luminosity 
of a cooling NS by Erot, we get the following expression for the cooling efficiency 



Lcooi ^ I 0-15 {Pqa/tq)'^ if core neutron - pairing and tq > 0.3, 

Ej.gf- [ 2.5 Po.i/Tg if no core neutron — pairing and r > 1. 



6^ 



Prom expressions ( |67| ) and (69) one can see that the pulsars with ages > 10® yr are most likely 



candidates for those with luminosities dominated by the PC heating. 



4.2. Numerical calculation of polar cap heating luminosity 



The returning positrons accelerate through the same potential drop as the primary electrons 
and their energy heats the PC surface. We can calculate the luminosity heating the PC by using 
expression (62) and substituting the corresponding quantities from the numerical calculation. 
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Figure 6 shows the dependence of on ^ (from the integrand of Eqn. |£2[)' which reflects the 

distribution of heating across the PC. Although the distribution of /9+(zo,i^) has maxima at = 
and 1^ = 1 (see Figure 3), 7inax(C) has a strong maximum at = and decreases monotonicahy 
with ^ (see Figure 4 and 5 of HM98), causing Lj^ to decrease with ^. 

The total positron heating luminosity, scaled with the spin-down luminosity as a function 
of characteristic pulsar age, P/2P is shown in Figure 7. The numerically computed Lj^/Erot 
increases nearly linearly with r in the same way as the analytic estimate in the unsaturated regime 
(first expression of Eqn. ||6^). Lj^/Erot also increases approximately linearly with period for a 
constant age. The PC heating luminosity is thus a negligible fraction of the spin-down luminosity 
in young pulsars with ages less than r ~ 10^ yr, but becomes a significant source of heating in 
older pulsars. The numerical values are a factor of ~ 6 lower than the estimated values of /E^ot 



from the first expression of Eqn. (37) for the normal pulsars. For the ms pulsars, the analytic 
values are somewhat higher than the numerical values, but by less than a factor of ten. The reason 
that the analytic estimate, which simply assumes that the returning positron flux is half of the 
charge deficit at the PFF, is closer to the numerical values for normal-period pulsars may be that 
the screening scale heights are small and the returning positron fraction is not very sensitive to 
the structure of the electric field in the screening region. On the other hand the screening scale 
heights for the ms pulsars are about an order of magnitude larger, refiecting the larger (see Figure 
4) growth scale of the pair source function. Consequently, the returning positron fraction for ms 
pulsars is more sensitive to the structure of the electric field in the screening region. 



Although the analytic estimate of Eqn. (|67| ) seems to be a very good estimate of the returning 
positron luminosity, we caution that it cannot be extrapolated to ages much beyond the numerical 
results, i.e. above r ~ 10^ yr for normal pulsars and above r ~ 10^ yr for ms pulsars. This is 
because the assumption of complete screening breaks down at age-period combinations where the 
pulsar cannot produce enough pairs. When the pair density grows too slowly and cannot reach 
a high enough level for complete screening, then the returning positron fraction drops below the 
charge deficit prediction. Beyond the upper end of each constant period line in Figure 7, numerical 
solutions of returning positron fraction and screening scale length do not exist. In fact, complete 
CR screening terminates on the constant period line at about the same point where PFFs are no 
longer produced. Beyond this point, the returning positron fraction and the heating luminosity 
will drop sharply. However, as we will discuss in our next paper, ICS pair fronts are produced by 
older pulsars, even when the ICS screening from these pairs is incomplete. This will result in a 
PC heating component due to returning ICS positrons from partial screening which can be quite 
significant for ms pulsars. 

The results in Figure 7 have assumed a constant value of inclination angle, x = 0-5. The 
expected variation of /+ oc cos^x (where p = 5/7 or 1/2) will cause /+ (= Lj^/Erot) to gradually 
decrease with increasing x until x approaches 7r/2. At this point, the second term in Eqn. (|5|) 
for the electrostatic potential (the relativistic equivalent of A83 solution) becomes important. For 
X = vr/2, our solution for the PC heating rate will be comparable to that of Eqn. (|6^). 
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X-ray emission at energies 0.1-2 keV has been detected from several dozen pulsars by 
ROSAT (Becker &: Trumper 1997). There is a rough empirical correlation of X-ray luminosity 
with spin-down luminosity, giving ^ 10^'^ Eroti so that the observed level of X-ray emission in 
pulsars would be a constant line in Figure 7 at 10"^. The maximum calculated values of / Ej-ot 
do reach the observed level, indicating the emission from PC heating in normal pulsars with 
r > 10^ yr and in ms pulsars is detectable. 

In Table 1, we compare our computed values of flux and surface temperature from PC heating 
with measured values of some older pulsars in which hot thermal components have been detected. 
Our values are those at the NS surface and have not been corrected for gravitational redshift effects 
- i.e. L+ and would be about 40% lower (a factor of o^) and T+ would be about 20% lower 
(a factor of a) for observers at infinity. Greiveldinger et al.(1996) have fit three-component (two 
thermal and one power-law) spectra to combined ROSAT and ASCA data from PSR0656-I-14 and 
PSR1055-52. These pulsars are young enough to have expected cooling as well as heating thermal 
components, so it is important to separate the emission from these two components in order to 
make comparisons with PC heating models. PSR0656-f 14, however, has an inferred surface dipole 
field of 9 X lO""^^ G and so is not in the regime of CR pair fronts near the NS surface that wc have 
treated in this paper. It may still have a CR pair front at higher altitude if the ICS pair front 
is unstable, a situation we will address in our next paper. The fits for the heated area A of the 
hot components in both cases are much smaller than the NS surface area, and are thus consistent 
with emission from heated PCs. Our predicted values of PC temperature from returning positrons 
agree fairly well with the measured values for these pulsars, although our predicted fluxes are 
significantly lower than those observed. Because these middle-aged pulsars are still dominated by 
cooling components and have power law components at high energies, extracting the relatively 
small heating component is very difficult. 

Wang & Halpern (1997) have fit single-component blackbody spectra to ASCA data of 
PSR1929-I-10. This pulsar is too old to have detectable cooling component and the one-component 
fits of ASCA data indeed indicate emission from heated PC whose heated area is much smaller 
than even the standard PC area. We have thus assumed the measured value of area A in 
computing the theoretical PC temperatures, r+ = (L^/aA)^/'^, of these pulsars, instead of the 
canonical PC area. Ape = ttR{VIR/ c). Our computed temperatures are in agreement with the 
measured values, within the uncertainties. Long period (older) pulsars have smaller PCs and thus 
the returning positron luminosity will heat a smaller area and we therefore would predict that 
pulsars with longer periods should have higher PC temperatures. PSR1929-I-10 does in fact have 
higher measured temperature and smaller heated area. Although Geminga is expected to have 
both cooling and heating components, it is more difficult to compare our theoretical results for 
positron heating with the observed values from a single-component thermal fit. The fact that 
our predicted flux is somewhat lower than the measured value is consistent with an additional 
cooling component. We have used the canonical PC area to compute our predicted temperature 
for Geminga, so we would expect it to be higher than a measured temperature which includes 
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cooling from the whole NS surface. 

5. SUMMARY AND CONCLUSIONS 

In this paper we derived some practical formulae for the positron fluxes returning from the 
upper PFF in rotation-powered pulsars in cases where E^^ screening is produced by pairs from CR. 
We presented the expected theoretical values for the PC X-ray luminosities due to the heating by 
precipitating positrons. The calculated efficiencies of PC heating explicitly depend on the pulsar 
spin-down age and spin period, and can be used in the analyses of thermal X-ray fluxes from 
pulsars. Our numerical calculation of returning positron flux and PC heating efficiencies show that 
the dependence on r and P is the same as in analytic formulae, the actual values of the numerical 
quantities are smaller than the analytic estimates by a factor of 2-3. In summary, we have reached 
the following conclusions. 

• The heating of the PC by returning positrons is possible if stable pair fronts develop, and E^^ 
is not assumed to be fully screened at the onset of the PFF. In most cases we discuss in this 
paper the fraction of returning positrons is not affected by the details of our modeling of the 
screening of E^y and is mostly determined by the distribution of pairs beyond the PFF. 

• In contrast to the results obtained earlier by A81, we find that the returning positrons may 
significantly contribute to the PC heating and therefore to the observed pulsar thermal 
X-ray fluxes, especially for older pulsars. 

• Our theoretical model can be thoroughly tested against observations as soon as more data 
on the pulsed thermal X-ray fluxes from middle-aged and old pulsars become available. We 
anticipate that the pulsar X-ray light curve would have larger pulse fraction in the case 
of NS thermal emission produced by the PC heating than in the case of thermal emission 
during NS photon-era cooling. 

• We predict that long-period pulsars will have higher surface temperatures from PC heating 
than those of short-period pulsars. 

• Our calculations indicate that the PC heating increases toward the magnetic axis, which 
could manifest itself in a fine structure of the X-ray light curves and heated areas smaller 
than the PC area. 

We caution that our analytic expression in Eqn. (|6^) for /+ is not applicable for normal 
pulsars with r ~ 10'' yr or ms pulsars with r ~ 10^ yr. The reason is that the analytic formulae are 
based on constant magnetic field and do not take into account the cessation of pair formation by 
CR photons beyond a certain age (dependent on period). Thus, our plots (see Figure 7) for pulsars 
with 5 and 10 ms periods cannot be compared with the observational data for actual ms pulsars. 
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most of which have ages r ~ 10^ yr where CR photons cannot produce pair fronts. Although our 
analytic formula for would give very large heating luminosities for these ms pulsars (beyond 
the region of validity of the formula), the actual heating efficiencies will be much lower. However, 
as we will address in our next study, ms pulsars are capable of producing pairs via non-resonant 
ICS and the resulting heating luminosities will be detectable. Furthermore, the results presented 
in this paper are not applicable to pulsars with surface magnetic fields above about 4 x 10^^ G, 
which are capable of ICS screening. The possibility of ICS pairs from those pulsars which are 
beyond their CR deathlines, as well as ICS screening and PC heating in higher-field pulsar, will 
be investigated in our next paper. 

We thank the anonymous referee for careful reading of the manuscript and helpful comments. 
This work was supported by the NASA Astrophysics Theory Program. 



- 25 - 



APPENDIX 

SOLUTIONS WITH EXPLICIT BOUNDARY CONDITIONS AT THE UPPER 

BOUNDARY 

Here we present the solutions to Eqn. (1) subject to the expHcit boundary conditions 
(especially at the upper boundary) different from those discussed in the main text. These solutions 
illustrate how the onset of the upper boundary affects the entire distribution of electrostatic 
potential and accelerating electric field above the PC. It is also instructive to compare these 
solutions with those we used in our calculations (see the main text for details). Throughout this 
Section we denote the dimensionless height of the upper boundary by Zc {= rjc — I). We factorize 
the expressions for the electrostatic potential and field by $0 = {^R/c)BqR and Eq = ^q/R, 
respectively. Finally, the solutions we present here imply the space- charge limitation of current. 

The solution for radial function J^i [see Eqns (^)-(|7|)] for which the following Dirichlet 
boundary conditions 

= 0) = ^{z = Zc) = 



are satisfied, is 



Z - Zc 



sinh^(7^z/2) 
■sinh^(7j2:c/2) 



smh[-ii{zc - z)] + sinh(7jz) - sinh(7iZc 
1 - cosh(7j2;c) 



(^1) 



Let us also present some asymptotic expressions. For the case where z < Zc « r^^R we get 



2 



-E^^z (1-24 



1 



Kcos X + -0oC-f^(l)'^(l) sin X cos > 



K cos X + 2^'^^^("'^-''^('''-' ^™ ^ ' 



(^2) 
(A3) 



For z « rpc/R < z^ we arrive at formulae (|20|), (|2lD . Finally, for z ~ rpc/R we get 

z 



2 ° c /(I) 



:i-a--8E#7|Texp[-7.(^e-.)] 



cos x+ 



1 



^0^^(1)5(1) 



c(i -e)-- lef: exp[-7.(^c - z)] 



sm X cos ( 



E. 



i|-"2^°~7(T)^^"^ ^ 



1 



K COS X + ^^oC-f^(l)'^(l) sin X cos < 



Solution satisfying the Neumann boundary conditions -Ey (z = 0) = (z = Zc 
from the following radial function 



(A4) 



(A5) 

can be derived 



= [lime - 2/7^) cosh7(zc - z) + (2?7c - ??) sinh7(zc - z) + (2 - 77) sinh(7z)- 



7(r/ - 2/7^) cosh(7z) - (7r?c - 2/7) cosh(7Zc) - (2r?c - 1) sinh(7Zc) + 7 - 2/7+ 
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(7^r?c - l)^;sinh(7Zc) + 'jZcZ cosh{jZc)]/[yricSmh{'yZc 
For Zc << 0.1(rpc/i?)^ the solution can be approximated by 



cosh(7Zc) + !]• 



^11 



3 ^ Zc 
-$o — 
2 1-e 



1 



1 



-z 1 



2z 

3Zc 

z 

Zr- 



K COS X + 2^oCH (1)6(1) sin x cos < 



K cos X + 2 ^0'^-f^(l)^(l) sin x cos < 



(>16) 



(A7) 



(A8) 



For z < 0.3rpc/-R < Zc we again arrive at formulae (| 
solution reduces to 

" nR z z 



El 



(J) = — (J)^ 

4 ° c /(l)z, 

9 f^i? 1 z 

~~2 °~7(i)i: 



1 



1 



2z 

3zc 
z 

Zr 



whereas for Zc » 0.1(rpc/ii)^ the 



(A9) 



K COS X + ~;(^oiH (1)6(1) sin x cos i 



1 



K COS X + ^^o^-f^(l)'^(l) sin X cos i 



(AlO) 



In the above equations the correction factors /, H, and 6 accounting for the gravitational redshift 
effect read (see also HM98) 

3 



f{x) 



H{x) 



X 

e 
e 

X 



In 1 



+ 



1 + 



+ 



1 



2x 

Ze/2x + K/2x^ 



(All) 



(A12) 
(AU) 



X- (1 - €/x)f(x) 

6(x) = d\n[H(x)e(x)]/dx, 

where x > 1, and 6(x) is the half-opening angle of the polar magnetic flux tube defined right after 
Eqn. (3). 



Note that in the case Zc « Tpc/R the imposing of the boundary condition ^(zc) = 
additionally supresses the electric field near the stellar surface [cf. Eqns (A3) and (A8)] compared 
to the case with the boundary condition E\\(zc) = 0. This is a clear illustration of the fact that 
^(zc) = condition forces E\\ to vanish well below Zc- This effect is especially pronounced when 
the upper boundary is very close to the stellar surface (e.g. when Zc « rpc/R). 

By comparing expressions (A2), (A3) with (A7)-(A10) one can see that the altitudinal offset 
between the vanishing of electrostatic potential and vanishing of parallel component of the electric 
field ranges from 0.3 to 0.5 Zc- This supports our conclusion (see §2.3) that simultaneous vanishing 
of both £^11 and Vy • E^^ hardly occurs within the PFF. Note also that Eqns (A5) and (A8) are the 
same as Eqns (A5) and (Al), respectively, presented in HM98. 

Recently, Dyks &: Rudak (2000) explored some useful approximations to -Ey that correspond 
to our general solution (A6). In particular, they arrive at the same expressions as (A8) and (19) 
[see their Eqns (9) and (12), respectively]. Here it is worth mentioning that in our previous paper 
(HM98) one of the approximate formulae, Eqn. (A3), is erroneous and should be replaced by e.g. 
formula (21) of this paper or by similar fitting formula given by Dyks & Rudak. Additional study 
aimed at the fitting of numerous combersome analytic expressions by simple compact formulae 
would be desirable. 
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Fig. 1. — Example of pair source function, integrated over energy, as a function of distance, x, 
above the PFF. The vertical axis measures the number of pairs produced by CR photons in each 
spatial bin, normalized per primary electron. 
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Fig. 2. — Self-consistent solution for the charge density, p, due to returning positrons which 
asymptotically approaches the charge deficit, Ap, needed to screen the electric field, E\\ , above the 
PFF, modeled as a declining exponential with screening scale height, A^i?. The pulsar parameters 
are the same as those of Figure 1. 
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Fig. 3. — Solutions for the returning positron density, p^/pcj, normalized to the GJ density and 
the screening scale height, A^i?, as a function of magnetic colatitude, ^, which has been normalized 
to the PC half-angle. 
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Fig. 4. — Solutions for the screening scale height, A^, as a function of surface magnetic field 
strength in units of the critical field, B/Bcr, for different pulsar periods. 
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Fig. 5. — Solutions for the returning positron density, p^/pGJ, normalized to the GJ density as a 
function of surface magnetic field strength in units of the critical field, B/Bct, for different pulsar 
periods. 
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Fig. 6. — Example of the variation of PC heating luminosity, L_|_, as a function of magnetic 
colatitude, ^, which has been normalized to the PC half-angle. The normalization of the vertical 
axis is arbitrary. 
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Fig. 7. — PC heating luminosity, L_|_, normalized to the spin-down energy loss rate, Erot, as a 
function of the characteristic spin-down age, r = P/2P, for different pulsar periods, as labeled. 
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Heated Pulsar Polar Caps 





0656+14 


1055-52 


1929+10 


Geminga 


P(s) 


0.384 


0.197 


0.226 


0.237 


r (yr) 


1.1 X 10^ 


5.6 X 10^ 


3.2 X 10^ 


3.4 X 10^ 


$h (erg cm"^ s-^) 


2.4+^-^ X 10-12 


l.Ot^;^ X 10-13 


1.7 X 10-13 


4.78 X 10-12 


A (cm^) 
rh(K) 


2.5li:? X 10'' 
1.5l°| X 10^ 


l.ll^:^ X 10^ 
S.ltl^ X 10^ 


3 X 10^ 
5.14l°:|3 X 10^ 


5.6t°| X 10^ 


L+ (erg/s) 


4.6 X 10^° 


5 X 10^0 


4 X 1030 


6.8 X 1030 


(erg cm~^ s"-^) 


1.5 X 10-13 


2.3 X 10-1^ 


5.6 X 10-13 


2.3 X 10-12 


r+ (K) 


2.6 X 10^ 


2.3 X 10^ 


6.9 X 10^ 


2.5 X 10^ 



Table 1: Measured values of hot thermal component flux heated area A and temperature 
Th of PSR0656+14 and PSR1055-52 are from Greiveldinger et al.(1996) and of PSR1929+10 is 
from Wang & Halpern (1997). Measured values for Geminga are for the total thermal component 
(Halpern &; Wang 1997). In computing from L+, we have assumed a solid angle of Att. 



